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Abstract 

We review the spin radical pair mechanism which is a promising explanation of avian 
navigation. This mechanism is based on the dependence of product yields on 1) the 
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hyperfine interaction involving electron spins and neighboring nuclear spins and 2) the 
intensity and orientation of the geomagnetic held. This review describes the general 
scheme of chemical reactions involving radical pairs generated from singlet and triplet 
precursors; the spin dynamics of the radical pairs; and the magnetic held dependence of 
product yields caused by the radical pair mechanism. The main part of the review includes 
a description of the chemical compass in birds. We review: the general properties of the 
avian compass; the basic scheme of the radical pair mechanism; the reaction kinetics 
in cryptochrome; quantum coherence and entanglement in the avian compass; and the 
effects of noise. We believe that the quantum avian compass can play an important role 
in avian navigation and can also provide the foundation for a new generation of sensitive 
and selective magnetic-sensing nano-devices. 
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1 Radical Pair Mechanism 


1.1 Introduction 

A radical is an atom, molecule or ion that has unpaired valence electrons. Radicals and 
radical pairs often play a very important role as intermediates in thermal, radiation, and 
photochemical reactions [1J. The presence of unpaired electron spins in these systems 
allows one to influence and control these reactions using interactions between external 
magnetic fields and electron spins [2]. However, until 1970, most scientists believed that 
ordinary magnetic fields had no significant effect on chemical or biochemical reactions, 
since the magnetic energy of typical molecules, under ordinary magnetic fields, is much 
smaller than the thermal energy at room temperature and is much smaller than the acti¬ 
vation energies for those reactions HU2J. This situation changed significantly in the 1970’s 
after a series of experimental results were reported on magnetic held effects on chemical 
reactions [3H7]. Because of these experimental studies, a number of researchers have made 
an effort to theoretically explain the magnetic held effects on the chemical reactions mm- 
Thanks to these and the subsequent efforts, we are now able to explain systematically 
magnetic held effects in terms of the radical pair mechanism. The radical pair mecha¬ 
nism was then successfully applied to explain the chemically induced nuclear polarization 
and electron polarization, which were shown to be based on the spin dynamics of radical 
pairs (2j. 

According to the radical pair mechanism, an external magnetic held affects chemical 
reactions by alternating the electron spin state of a weakly coupled radical pair, which 
is produced as an intermediate. The basic scheme of chemical reactions through the 
radical pairs is shown in Fig. [0 Radical pairs are usually produced as short-lived inter- 
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Figure 1. Reaction scheme of radical pairs generated from singlet and triplet precursors. 
Singlet and triplet radical pairs are represented by 1 [... ] and 3 [... ], respectively. kso,P,E are 
the rates of reactions. (Taken from Ref. [2], with modifications.) 
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mediates through decomposition, electron transfer, or hydrogen transfer reactions from 
singlet or triplet excited states. These reaction precursors are called “S-precursors” or 
“T-precursors”. 

The generated radical pairs are surrounded by a solvent molecular cluster, called a “sol¬ 
vent cage”, and these pairs retain the spin multiplicity of their precursors. Initially, two 
radicals are close together, and are called the “close pair”. Sometimes, recombination 
reactions occur from S- or T-close pairs immediately after the formation of the radical 
pairs. Such reactions are called “primary recombinations”, and their products are called 
“cage products”. However, because of the Pauli exclusion principle, the T-close pairs re¬ 
quire enough energy to produce excited “cage products”, while the S-close pairs are able 
to produce the ground state of the “cage products”, due to the spin preservation during 
the chemical reactions. Consequently, the recombination reactions from T-close pairs oc¬ 
cur less frequently than those from S-close pairs. Usually, we can ignore reactions from 
T-close pairs. However, the singlet-triplet conversion is possible for close pairs involving 
heavy atom-centered radicals due to their spin-orbit interactions. But for close pairs in¬ 
volving only light atom-centered radicals, no spin conversion occurs between their singlet 
and triplet states. The S-T conversion mainly occurs during the second stage. In the 
second stage, as shown in Fig. Q] the two radicals begin to diffuse away from each other, 
forming a separate pair. When the two radicals are separated at a certain distance, the 
S-T conversion becomes possible through weak magnetic interactions of radicals including 
Zeeman and the hyperhne interactions, as will be explained in detail later. In the last 
stage, some of the separated radicals approach each other, forming close pairs again, and 
some continue to diffuse from each other, forming free radicals and producing “escape 
products” with or without the solvent molecules DEI 
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1.2 Spin Dynamics of the Radical Pairs 

Consider two weakly coupled radicals that form a radical pair. The spin dynamics of the 
radical pair is governed by a Hamiltonian (Hpp), which can be expressed as the sum of 
an exchange term (H ex ) and a magnetic ( H mag ) term [2j[8], 


Hpp H ex -\- H ma g , (1) 

where, 

H ex = -J(r)(2S 1 -S 2 + ^j, ( 2 ) 

a b 

Hmag = IJ>bB • (fi'lS'i + g 2 S 2 ) + ^ A\jS\ ■ In + A 2 kS 2 ■ l2k) ■ (3) 

i fc 

In Eq. |5i = |<Tj, where i = 1,2, and o % are the Pauli matrices, J(r ) is the value of the 
exchange integral between two unpaired electron spins (Si and S 2 ), which decreases with 
separation distance, r. In Eq. [3l the first two terms describe the Zeeman effects, and 
the last two terms are hyperhne interactions between the electron spins (Si, S 2 ) and the 
nuclear spins (In, I 2 k) in the radicals 1 and 2. Nuclear Zeeman effects are neglected since 
their magnitudes are much smaller than those of the electron Zeeman terms and hyperhne 
coupling terms. Also, gi and g 2 are the isotropic g-values of the two component radicals 
in the radical pair, respectively, and An and A 2 k are the isotropic hyperhne coupling 
constants in radicals 1 and 2, respectively, and the number of nuclei in radicals 1 and 2 
are a and b, respectively. 

The state of a radical pair can be represented by the product of the electron and nuclear 
states. The two unpaired electron spins generate the singlet (IS 1 )) and triplet (| T n ); n = 
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Figure 2. Dependence of a radical pair’s energy on the distance (r) between two components 
when B ^ 0 and J(r) is negative. (Taken from Ref. |T], with modifications.) 
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- 1 , 0 , 1 ) states which are expressed as: 


\S) = 

^(iw)- in»), 

(4) 

\To) = 

T(m)+ in))), 

(5) 

T + 1 ) = 

itt), 

( 6 ) 

T- 1 ) = 

in). 

(7) 


We use |t) or |t) to express the ^-component of a electron spin state, and we use | M) to 
express the ^-component of a nuclear spin state, i.e. I Z \M) = M\M). Since we ignore 
the interactions between nuclear spins, the total nuclear spin states can be expressed as, 
\Xn) = II i l^i) rifc where a and b are the number of nuclei in radicals 1 and 2, 

respectively. 

In the presence of an external magnetic field, we can calculate the energies of the singlet 
and the triplet radical pairs, 

E(S) = (S, X n\H rp \S,Xn) = J(r), (8) 

a b 

E(T n ) = (T n ,XN\H RP \T n ,XN) = - J(r ) + ngn B B + 77 A u Mj + ^ A 2 k M k ) , (9) 

i k 

where, n — - 1,0 and + 1 ; g = 9l + 92 . 

The r-dependent energy of a radical pair in the external magnetic field is schematically 
depicted in Fig. [2] in the range of a few molecular diameters. Since the exchange interac¬ 
tion (J(r)) decreases with the distance exponentially, J(r) can be neglected safely in the 
second stage (r > r c ) of radical pair mechanism mentioned before. Therefore, for further 
discussion, we will neglect the exchange interaction (J(r)). 



When r > r c , the conversion between | S) and \T n ) is governed by the following rules: 


^ a b 

(To, Xn\Hrp\S, xn) = 2 (^QUbB + A\ t M t — A 2 k M k )) = Tiujn, (10) 

i k 

{T±i,Xn'\Hrp\S, xn) = ^■[h i (hi + X ) “ M i( M * T !)]*> ( n ) 

where, Ag = g x - g 2 , \ Xn) = \ M u M k ), \Xn’) = \M',M' k ), M[ = Mj =f 1, M' k = M k , and 
u>n is the nuclear Rabi frequency. 

It is worth mentioning that in large external magnetic fields, due to the large energy 
gap (gfisB), only the S — T 0 conversion is taken into account, and the S — T± i conver¬ 
sions can be neglected m However, in small magnetic fields, all conversions between 
singlet and three triplet states of the radical pair must be considered, since the Zeeman 
terms are small compared to hyperfine interaction. 

By solving the Schrodinger equation (fl2lh Kaptein obtained the time evolution of the 
wave function ($(i)) of a radical pair during the S — To conversion [IT]. 

= H RP m, ( 12 ) 

$(t) = CsN(t)\S, xn) + Ct 0 n\Tq, Xn), (13) 

where, CsN(t ) = Cs(0) cosoj^t—iCT o (0) sin uj^t, Ct 0 n(T) = Cr o (0) cosoj^t—iCs(0) sin uj N t. 
For chemical reactions occurring from S-precursors, |C's r iv(0)| = 1 and |Ct o v(0)| = 0. In 
these cases, the wave function of a radical pair is, <h(t) = (cos cojsrt)\S, Xjv)+(sino;jvt)|2o, Xn) 

Ideally, the time evolution of the populations of singlet and triplet states are cos 2 uj^t 
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and sin 2 cj during the S — T 0 conversion, respectively. However, due to random-walk 
diffusion within the solvent cages, the radicals either separate from each other further, 
forming “escape radicals”, or re-encounter with each other, realizing a secondary recom¬ 
bination. Therefore, the solvent can affect the lifetime of the generated radical pair. The 
more viscous the solvent is, the longer the lifetime of the radical pair becomes. 

1.3 Magnetic Field Effects on Product Yields due to the Radical 
Pair Mechanism 

The magnetic field can have an influence on the products of chemical reactions through 
the radical pairs in solution [12] . Since the lifetime of a radical pair becomes longer in a 
more viscous solvent, the magnetic field effects, due the radical pair mechanism, become 
appreciable if the solvent is very viscous [3]. 

Eq. (ROD and Eq. (fTT| clearly show that the S — T conversion of a radical pair can 
be influenced by the external magnetic field and the hyperfine interaction, in the region 
where the exchange interaction (J(r)) can be neglected. Thus, the magnetic field effects 
on radical pairs in this region can be classified by two typical mechanisms: the A g mech¬ 
anism and the hyperfine coupling mechanism. The A g mechanism is appreciable when 
A g ^ 0 and An = A 2k = 0, while the hyperfine coupling mechanism is appreciable when 
A g = 0, and A u ^ 0 or A 2k ^ 0 [1][2]. 

For the A g mechanism, there is no S — T conversion in the absence of external mag¬ 
netic field. However, according to Eq. (ITU , in an external magnetic field, the S — To 
conversion can be induced with a frequency of a 9i j -b R / 2 h, but the S — T ±i conversion will 
still not occur. This means that, due to the A g mechanism, the rate of S — T conversion 
becomes larger with a stronger external magnetic field. In this situation, if the radical 
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pair in a chemical reaction is produced from an S-precursor, the cage product will de¬ 
crease as the magnetic field increases. Because a singlet radical pair is generated initially, 
the conversion rate will increase as the field increases. Also, it can be proved that the 
magnetically induced decrease in the yield of the cage product is proportional to B I pQ. 
On the other hand, the yield of products from the triplet radical pair, such as escape 
products, increases with B increasing. The opposite results occur when the radical pair 
is generated from a T-precursor. The yield of cage product increases, with B increasing 
at a rate proportional to B^. The yield of products from the triplet pair, such as escape 
products, decreases. 

For the hyperfine coupling mechanism at zero magnetic field, the S — T conversion can 
occur between the singlet state and all of the three sub-levels of the triplet state, through 
the hyperfine coupling terms in Eq. (TTOT) and Eq. (ITT]) . In the presence of an external 
magnetic field, the three triplet sub-levels are split by the Zeeman effect. So the S — T± i 
conversion will be forbidden due to the energy gap, but the S — T 0 conversion can occur 
with the same frequency, (E“ A u M % - '± 2 l A 2 k M k)/ 2 h, as the frequency without the magnetic 
field. As a result, the rate of the S — T conversion decreases as the external magnetic field 
increases. And this decrease becomes saturated at higher fields. Due to this mechanism, if 
the radical pair is produced from the S-precursor, the yield of cage products will increase, 
and the yield of products produced from the triplet pair will decreases, when a non-zero 
magnetic field is present. The magnetically reduced triplet yield is characterized by the 
field, T>i/ 2 , at a half saturation. Weller et al. [13] gave an expression of the Bi/ 2 value as, 
By, = \{B X + B 2 ), where B t = + X ))^ (! = 1 or 2 ) 3 = l, 2 , •••>« for / = 1 

and j — 1, 2,..., b for l — 2). The opposite occurs for a pair generated from a T-precursor. 
The yield of cage product decreases with B increasing at a rate proportional to B^. The 
yield of products produced from the triplet pair, such as escape products, increases with 
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B , and these magnetically induced changes should also become saturated at very high 
magnetic helds (B Bi/ 2 ). 

It is noteworthy that due to the A g mechanism and the hyperfine coupling mechanism, 
the external magnetic held has opposite effects on the product yields. The general case 
will be a mixture of these two mechanisms, meaning A g ^ 0, and A lt ^ 0 or A 2 k ^ 0. In 
this case, the magnetic field effects, due to the A g mechanism and the hyperfine coupling 
mechanism, appear simultaneously. Since the magnetic field effects due to the hyperfine 
coupling mechanism have a saturated region, the magnetic field effects due to the A g 
mechanism dominate at high fields, and the magnetic field effects due to the hyperfine 
coupling mechanism dominate at low fields. Therefore, the cage product yield has a max¬ 
imum (nadir), when the radical pair is generated from an S-precursor (T-precursor) as a 
function of increasing magnetic field. 

2 Chemical Compass in Birds 

2.1 Introduction 

Before the twentieth century, biology and physics rarely crossed paths. But after the 
development of quantum mechanics, Schrodinger introduced quantum physics into the 
realm of biology [2] • However, quantum mechanics normally applies at the microscopic 
regions of the order of nanometers, and living organisms belong to the macroscopic world. 
Usually, quantum collective (coherent) effects can only be observed in laboratories using a 
high vacuum and/or an ultra-low temperature. Despite these differences, the progress in 
quantum biology has been very rapid since the twentieth century [29 j l44 ] l45 ] l48]. In some 
sense, the quantum mechanics always plays a significant role in biology since every chem¬ 
ical reaction relies on quantum mechanics, and chemical reactions are the basic processes 
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in biological systems [T5HT9] . Recently, many articles have been published in the field 
of quantum biology, including articles on photosynthesis [2Ul - f39] . avian compass mm, 
and olfaction [561157] . 

The precise mechanism of avian navigation has been a mystery for centuries. Every 
year migrant birds navigate hundreds or even thousands of miles between their seasonal 
habitats. Scientists have been trying to explain this astonishing phenomenon for decades. 
To date, the primary biophysical process underlying the avian compass still remains un¬ 
known. Scientists found that many sources of directional information can contribute 
to birds’ navigation, such as landmarks or topographic features [5B1EH], stars fClOl 1(TT1 . 
the sun [6ol 1671 . and so on. Among these sources and hypotheses, magnetite as well as 
biochemical radical-pair reactions has been hypothesized to mediate sensitivity to Earth- 
strength magnetic fields through fundamentally different physical mechanisms [681169] , 
In the magnetite-based mechanism, magnetic fields exert mechanical forces Mm- Re¬ 
cently, Treiber et al. reported that the iron minerals (located in cells in the upper beak of 
pigeons, and possibly magnetite), are actually macrophages and have nothing to do with 
magnetic sensing [721173] . but the magnetite-based mechanism is still a reasonable hypoth¬ 
esis to explain the magnetoreception in birds. In the radical-pair mechanism, magnetic 
fields alter the dynamics of transitions between spin states, after the creation of a radical 
pair through a light-induced electron transfer. These transitions in turn affect reaction 
rates and products (SUSHUSB] . 

2.2 The Properties of the Avian Compass 

There are also two important properties of the avian compass that can be explained by 
radical pair mechanism very well. 
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First of all, the avian compass is an inclination compass rather than a polarity com¬ 
pass (Z3HZZI. The geomagnetic field lines leave the earth at the magnetic pole near the 
geographic south pole, curve around the earth and re-enter to the earth at the magnetic 
pole near the geographic north pole. The inclination, defined as the angle between the 
field lines and the horizontal plane, is 90° at the magnetic poles and 0° at the magnetic 
equator; at the latter, the field lines run parallel to the earth’s surface. In the northern 
hemisphere, field lines have a downward component with a positive inclination; in the 
southern hemisphere they have an upward component. Therefore, the geomagnetic field 
provides a reliable, omnipresent source of navigational information for birds capable of 
sensing the field. The functional mode of the avian magnetic compass is different from 
that of our standard compasses. The avian magnetic compass was found to be an inclina¬ 
tion compass based on the inclination of the field lines instead of their polarity m Birds 
can only perceive the axial course of the field lines and they must interpret the inclination 
of the field lines with respect to up and down to derive unambiguous directional informa¬ 
tion. The avian magnetic compass does not distinguish between magnetic “north” and 
“south” as indicated by polarity. 

A second important property of the avian compass is its light dependence. This property 
has been supported by a series of behavioral experiments [T3UT81 - I84] . The avian compass 
orientation is dependent on the wavelength of the ambient light. While European robins 
and Australian silvereyes demonstrated good orientation in blue light and green light, 
they were disoriented in red light j411 ITS] [79]. Also, European robins cannot orient in 
yellow-orange light, although there is a significant overlap between the green light, where 
orientation was excellent, and the yellow light without orientation [ST]. These experimen¬ 
tal findings put constraints on any mechanism designed to explain magnetoreception in 
birds and provide very strong support for the radical pair mechanism [85] [86]. 
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2.3 Basic Scheme of the Radical Pair Mechanism 


The actual reaction scheme of the magnetoreception can be very complex in birds. It 
involves molecules in the eyes absorbing the light, the chemical reaction of the radical 
pair in these molecules, and the avian neurons to react to the signal. However, we can 
simplify the scheme into three stages. 

In the first stage, the photons with enough energy activate a certain type of molecules 
located in the birds’ eye, inducing an electron transfer reaction and generating the radical 
pairs in their excited singlet states. After the pair is generated, under the influence of 
the external magnetic field (the geomagnetic field) and the internal magnetic field (the 
hyperfine coupling effect), the state of the pair can remain a singlet state or become a 
triplet state. A different inclination angles, associated with the external magnetic field, 
can induce different ratios of the singlet and triplet states. In the last stage, the molecules 
in different states will generate different chemical products which can induce a detectable 
signal that the birds can use to recognize the direction they need to go 0. 

The predictions of the radical-pair mechanism are consistent with behavioral findings: 
the reaction yield is independent of the polarity of the magnetic field [751188] , reception of 
magnetic information takes place in the eye j89j[90], reception is strongly affected by the 
ambient light conditions [THflSTL 183] . and reception is consistent with the postulated role 
of ocular photoreceptors in creating magnetosensitive radical pairs [85j . 

Although there are not any observations of specific biological radical pairs satisfying all 
constraints, the current discussion of cryptochrome la (Cryla) makes it a promising can¬ 
didate for magnetosensitive radical pair [ 87llTnTl96] . Also, the electron transfer path has 
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been studied [971 [98] . However, the neural path is still unknown, although a vision-based 
hypothesis was put forward by Ritz et al. m- 

2.4 Reaction Kinetics in Cryptochrome 

Using a carotenoid-porphyrin-fullcrene model system [99], Hore and his colleagues demon¬ 
strated that, under a static magnetic field, as weak as the geomagnetic held, a chemical 
reaction can act as a magnetic compass, by producing the detectable changes in the 
chemical product yield. Their experiment has provided a proof-of-principle of a chemical 
compass. Along with other experiments [75 |[7HH82] , this provides a significant support to 
the radical-pair hypothesis of the avian magnetic compass. 

However, the candidate molecules with certain biophysical characteristics must exist in 
the eyes of the migratory birds, so that the radical pair mechanism can play a role in 
the magnetoreception. Cryptochromes, a class of photoreceptor proteins located in the 
retina and absorbing blue-green part of the spectrum which is the functional range of 
the magnetic compass [100 ] . were proposed as the host molecules for the crucial radical 
pair cofactors that putatively act as primary magnetoreceptors. Also, the experiments 
have demonstrated that the vertebrate cryptochromes, under the UV/visible transient 
absorption and the electron paramagnetic resonance spectroscopy, form long-lived radi¬ 
cal pairs which involve a flavin radical and a radical derived from a redox-active amino 
acid [101. 102] , This demonstrates that the cryptochrome harbors are the type of radi¬ 
cal pair needed for the action of the magnetic compass. Scientists also have verified that 
cryptochromes exist in the eyes of the migratory birds [9111100111011110311104] . Furthermore, 
a distinct part of the forebrain, which primarily processes input from the eyes, is highly 
active at night in the night-migratory birds [10511109] . All of these findings are consistent 
with the hypothesis that the cryptochromes can serve as the primary magnetoreceptors. 


16 










The process of cryptochrome photoactivation has been discussed, and several reaction 
schemes have been proposed [[971 11031 fll0fll20j . Cryptochrome contains a blue-light- 
absorbing chromophore, flavin adenine dinucleotide (FAD). This FAD cofactor is reduced 
via a series of light-induced electron transfers, from a chain of three tryptophans (Trp) 
that bridge the space between FAD and the protein surface [ 119j . 

According to Hore and his colleagues [98], due to the light-induced electron transfer, a 
radical pair [FAD*~TrpH* + ] (RP1) is formed in the cryptochrome (Cry-1 from the plant 
Arabidopsis thaliana, AtCry), followed by a second radical pair [FADH'Trp*] (RP2). The 
formation of RP2 is more complex. The tryptophanyl radical may deprotonate either 
fully or partially, while the protonation of FAD*” produces a the neutral FADH*. In this 
scheme, the RP1 interconverts coherently between singlet and triplet states, under the 
influence of the magnetic interactions internal to the radicals (electron-nuclear hyperfine 
couplings), and under the Zeeman interactions with the external magnetic field. Only 
the singlet state of RP1 can revert to the ground state [FAD+TrpH] by the electron-hole 
recombination, and the corresponding reaction of the triplet state being spin-forbidden. 
Simultaneously, one of the constituents of the RP1 changes its protonation state to give 
RP2, a process that is not spin-related, and in which the singlet and triplet are produced 
at equal rates. 

2.5 Applications Inspired by the Radical Pair Mechanism 

Despite all of the theoretical arguments, the ultimate goal of studying the mechanisms 
of bird navigation is to learn from nature and to design highly effective devices that can 
mimic biological systems in order to detect weak magnetic fields, and to use the geo¬ 
magnetic field to navigate. Previous literature has shown that the anisotropic hyperfine 
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coupling plays a crucial role in the magnetic field sensitivity of the avian compass and 
has provided routes towards the design of biologically inspired magnetic compass sen¬ 
sors [12F T221I124H131] . Scientists have exploited many practical methods to realize a 
device based on the avian chemical compass. Recently, several intriguing models were 
reported, and significant experimental/simulative results were obtained. 

One of the models is a synthetic donor-bridge-acceptor compass, which is a triad composed 
of linked carotenoid (C), porphyrin (P) and fullerene (F) groups |99lll22[fl23U132j . Such 
a triad molecule is the first known molecule that has been experimentally demonstrated 
to be sensitive to the geomagnetic field, although it works at low temperature (193K). 
In such a molecule, both radicals are immobilized. The anisotropic magnetic interac¬ 
tions are preserved because the hyperfine interactions in the fullerence radical are very 
small, and the anisotropic interaction in the carotenoid radical dominates the anisotropy 
of magnetic field effect. The other conditions are that the radical pair is generated in the 
singlet state with different decay rates for singlet and triplet states (fc-r < ks ), and the 
rate of the spin-lattice relaxation is comparable to the combination rate. Based on all 
such conditions, the experiments showed a promising result under a magnetic field of the 
order of ~mT and produced the surprising result that the lifetime of the radical pair is 
significantly affected by magnetic fields as weak as 50 /iT, which is about the magnitude 
of the geomagnetic field [99] . Yet, the fact is that currently there is not a biomimetic or 
synthetic chemical compass that functions at room temperature. 

Another model is one utilizing magnetic nanostructures to design a chemical compass 
inspired by theoretical studies [ 133] . According to the theoretical and experimental stud¬ 
ies, the internal anisotropic magnetic field is the key factor for product yields of the 
photochemical reactions to be sensitive to the angle with respect to the external mag- 
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netic field [Ml 03] IS]. In some molecules, the hyperfine coupling provides the internal 
anisotropic magnetic held. However, when designing a real device, the hyperfine coupling 
can be replaced by a local strong gradient magnetic held, which can be created in the 
vicinity of a hard ferromagnetic nanostructure by applying a spatially uniform bias held 
that cancels the held of the nanostructure in a small region [13311134] , In simulations, 
such a device exhibits signihcant directional sensitivity [133]. 

Other studies have achieved to map nanomagnetic held using a radical pair reaction |130j 
and found a method to improve the magnetic sensitivity of chemical magnetometers [131 j . 

2.6 Quantum Coherence and Entanglement in Avian Compass 

Biological mechanisms, such as magnetoreception in birds, are powered by the chemical 
machinery, which consists of complex molecules structured at the nanoscale and sub¬ 
nanoscale. The dynamics of the chemical machinery at such small scales is regulated 
by the laws of quantum mechanics. As one of the main features of complex quantum 
systems, the quantum coherence between different sub-systems, can ultimately lead to 
many collective phenomena, such as the quantum entanglement. 

2.6.1 Effects of Noise 

2.6.1.1 Thermal Noise 

Quantum effects, such as coherence and entanglement, are easily destroyed by interactions 
with the environment. There is a rule of thumb, the fcgT argument, stating that when the 
energy of the interaction is smaller than room temperature, then quantum coherent effects 
cannot persist. Intuitively, because the temperature of the biological systems is around 
300 K, which is hot, thermal fluctuations will destroy the quantum coherent effects. How¬ 
ever, the ksT argument can break down when dealing with living systems like birds. The 
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thermal argument is true only for the equilibrium state, which the system approaches for 
relatively long times. However, if the equilibrium is approached rapidly, coherent quantum 
effects can survive (at equilibrium) even at room temperatures [20l - f26lH35l[136j . Another 
possibility is that the quantum coherent effects, including the quantum entanglement, can 
survive as the system approaches the equilibrium with the surrounding environment. In 
another word, the quantum effects are used before the system has time to equilibrate with 
the environment. The latter situation is believed to be what happens in the radical pair 
mechanism of the avian compass [ 137] . 

Since the geomagnetic field is very weak 0.5 Gauss), the effects of such a weak mag¬ 
netic field, in biological systems, can easily be masked by the thermal fluctuations. To 
take advantage of the weak geomagnetic field, one effective method is to use the short¬ 
lived, specialized photochemical reactions. Then the thermal fluctuations do not have 
enough time to effectively mask the effects of the weak magnetic fields. The radical pair 
mechanism is believed to rely on this process [92] . 

If the useful signal carries more energy than the thermal noise, this signal can be dis¬ 
tinguished from the noise. But since the energy differences, due to the geomagnetic fields, 
between the singlet and the triplet states are generally small, and are even smaller in liv¬ 
ing organisms, the chemical reaction of radical pairs is a subject to a significant thermal 
noise [138] . However, due to the short lifetime of the radical pair reactions, and due to the 
fast conversion rate between singlet and triplet states, the photochemical reactions can 
produce a bird-detectable signal through the yields of singlet and triplet states, in spite 
of the thermal noise. Weaver et al. also discussed the detection limit in the radical pair 
mechanism, and claimed that a chemical compass in birds is feasible with the chemical 
reactions with certain rate constants 
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2.6.1.2 Dephasing Noise 


Other than thermal noise, dephasing noise can only affect the coherences between singlets 
and triplets. Once the radicals come close to each other, dipolar and exchange interaction 
can no longer be neglected. For instance, if the exchange interaction is considered, the 
electron singlet and triplet states will be separated in energy, which can destroy the con¬ 
version between the singlet and the triplet [ 140j . However, since the molecule motion can 
be quickly distributed over the environment owing to Brownian molecule motion, the mi¬ 
croscopic causes of the path taken are quickly lost. In addition, the involved energy scales 
are much smaller when compared with the energy released owing to recombination [140j . 
Therefore, the dephasing noise can hardly harm the coherence since the molecular mo¬ 
tions that cause spin-decoherence by modulating hyperhne interactions have either low 
amplitude or high frequency or both. 

There are many papers that discuss the effects of dephasing noise on the radical pair 
mechanism via Lindblad operators [125 :,lT40, il4l| . Intuitively, one might expect that de¬ 
phasing noise is unfavorable for the function of the chemical compass. However, these 
studies show that the effect of noise depends on the model. In Gauger’s model, the 
compass mechanism is almost immune to pure phase noise [141] , Cai argued that cor¬ 
related dephasing noise could even enhance the chemical compass in their model [125 ], 
Both models give us positive, or at least non-negative, perspectives concerning dephasing 
noise. 
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2.6.2 Quantum Coherence 


As we know quantum coherence describes a particle’s ability to exist simultaneously in 
several distinct states, such as position, energy, or spin. Quantum coherence plays a crucial 
role in the energy transport in photosynthetic complexes [24ll27LI29Lll42hll46lll48j . There is 
indirect evidence indicating that the quantum coherence exists in bird navigation [411147] . 
The proposed radical pair mechanism for avian compass postulates that absorption of 
blue or green light creates a radical pair of electrons. The radical pair undergoes coherent 
quantum oscillations between entangled singlet and triplet states, at a rate depending on 
the external magnetic fields. Finally, the pairs in the singlet and the triplet states will 
lead to different chemical products, providing orientation information for birds. So, the 
quantum coherence plays an important role in this mechanism. 

2.6.3 Quantum Entanglement 

Quantum entanglement is the quantum mechanical property that allows two or more 
particles to be correlated differently and stronger than the same classically correlated 
particles. This also means that while the global state of the system is perfectly known, 
the local state is fully mixed. The entanglement creates the non-local correlations, and 
the non-thermal excitations [137111470 l'49j . 

Quantum entanglement plays a key role in the radical pair mechanism of the avian com¬ 
pass. A single electron photo-excitation and a subsequent electron translocation leads to 
an entangled state, which provides the necessary spin correlations. Along with the effects 
of the magnetic fields, this causes the electrons to oscillate between singlet and triplet 
states. 

Many scientists have conducted researches on the role of entanglement in chemical corn- 
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pass fmi 1150111711 . Gauger et al shown that even though entanglement falls off at a 
faster rate than the decay of population from the excited state under dephasing noise, 
the superposition and entanglement are sustained for at least tens of microseconds, ex¬ 
ceeding the durations achieved in the best comparable man-made molecular system m- 
However, other researchers focus on the effects of the entanglement of the initial states of 
the radical pairs [15011171] . These studies have shown that the sensitivity of the chemical 
compass depends on the initial state of the radical pair. Briegel and his colleagues have 
randomly chosen 5000 different initial states from the set of separable states and singlet 
initial state and calculated the maximal achievable magnetic-field sensitivity for every 
value of B. The results show that entanglement is indeed helpful, and it is specifically 
entanglement rather than mere quantum coherence cm 

“Negativity”, originated from the observation due to Peres [151] . is a computable measure 
of the entanglement. It is based on the trace norm of the partial transpose, p TjX , of the 
bipartite mixed state density matrix, p |152]ll53j . Essentially, negativity measures the 
degree to which p TA fails to be positive, and therefore it can be regarded as a quantitative 
version of Peres criterion of separability [15114153] . From the definition of negativity, we 
can construct a quantity to measure the entanglement: N(p) = — where ||p Tj4 ||i 

is the trace norm of the partial transpose of the system’s density matrix [15311154] , This 
measure of entanglement, N(p), corresponds to the absolute value of the sum of negative 
eigenvalues of p TA [1 5511156] . and it vanishes for unentangled states. Moreover, in the case 
of an entangled state of two qubits, negativity is defined as two times the absolute value 
of the negative eigenvalue of the partial transpose of a state [157] . 

Denote p as a generic state of a bipartite system with the finite-dimensional Hilbert 
space, shared by two parties, A and B. p TA , the partial transpose of the density 
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matrix, p, with respect to subsystem A, will be defined as: 


(iA,jB I P Ta | rn A ,n B ) = (m A ,j B I P I i A ,n B ), (14) 

where | i A ,j B ) =| i) A ® \ J])b G 'H A ®T-L b is a fixed but arbitrary orthonormal product 
basis [153] • The trace norm of p TA is: 

|| p TA || 1 = try/ p Ta] P Ta , (15) 

which is equal to the sum of the absolute values of the eigenvalues of p TA , since p T - 4 is 
hermitian mu. Since the eigenvalues of the density matrix, p, are positive, the trace 
norm of p is: || p ||i= trp = 1. Thus, the partial transpose, p TA , also satisfies the 
condition: tr[p] = 1. But since it may have negative eigenvalues, e* < 0, its trace norm is: 
|| P Ta ||i= 1 + 2 | JAe* 1= l + 2A^(p) [153] . Therefore, the negativity (N(p) = \ Yhi e i\) can 
also be defined as the sum of the negative eigenvalues, +, of the density matrix partial 
transpose, p T ^, measuring by how much p TA fails to be positive definite [ 1531158] . So we 
have an alternative way to calculate the negativity. It can be written as: 

N (P) =1 ei l = ' Aj 2 Aj ’ 

* 3 

where e* are the negative eigenvalues of p Tj4 , and A j are the eigenvalues of p T + In the 
following, we will show some simple examples to illustrate negativity. 

Assuming the state \ijj) can be separated into two sub-states, | <f> A ) and \4>b), indicating that 
this state is not an entangled state in which \(p A ) = ^=(| t)+| I)) and \(p B ) = -j=(\ t) + | I)), 
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the expression of \ip) will be: 


— I 4>a) <e> | 4>b) — 77(1 T)a + 14-)a) <h> (| t)s + I 1)b)- 


(17) 


As a result, the density matrix of | -0), Pp, will be: 


Pp — 


|- 0 > (-01 

/ 


1 

4 


1111 

1111 

1111 

1111 


\ 


(18) 


Therefore, the partial transpose of the density matrix with respect to subsystem B is, 


Pp B 


1 

4 


/ 1 1 1 1 1 


1111 

1111 

1111 


V 


(19) 




Through the calculation, we can find that the four eigenvalues of p^ A are {1, 0, 0, 0}. Ac¬ 
cording to the formula of negativity above, N(p) = j ^4^, we find that the negativity 
is 0. This example shows us that the negativity of a unentangled state is zero. 


Next, we will consider an example of an entangled state | ijj) = ^(Iti) - lit)) = -^(0 11 0) 5 
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In this case, the density matrix of the state | -0), p^, is 


/ 0 0 0 0 \ 


1 01-10 


( 20 ) 



0-110 


y o o o o 


Therefore, the partial transpose of this density matrix with respect to subsystem B is 


/ 0 0 0 -1 \ 


T 1 0 10 0 
p? = 2 

z 0 0 10 


( 21 ) 


y -i o o o 


We can find that the eigenvalues of are {—0.5, 0.5, 0.5, 0.5}. Therefore, the negativity 


of the state 10) is 0.5. After taking twice the negativity, we find that the entanglement 
of the state |0) is 1, which is the maximum entanglement. 

The above examples show that negativity can serve as a metric of the entanglement. 
Therefore, in the following calculations we will use the negativity as the measurement of 
entanglement. However, it is worth to mention that there is another measurement of the 
entanglement named concurrence [159] , For a two-qubit system, the two measurements, 
negativity and concurrence, are equivalent. 

3 Calculations and Results 

In this section, we will introduce some of our work on the topic of the radical pair mech¬ 
anism [ 16011161] , We investigate two aspects of the radical pair mechanism: the yields of 
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Figure 3. Angular dependence of the triplet yields. The triplet yields are symmetric around 
90°. 

defined signal states and the entanglement of the states. 

3.1 Basic One-Stage Scheme 

Following Ref. [3T| , we include only the Zeeman interaction and the hyperhne interaction 
in the Hamiltonian of the system: 

2 

H = g f ji B Y,Si-(B + A i -fy. ( 22 ) 

2=1 

In Eq. 1)22]) . the first term is the Zeeman interaction and the second term is the hy¬ 
perhne interaction. (We assume that each electron is coupled to a single nucleus.) /, is 
the nuclear spin operator; S t is the electron spin operator, i.e., S = <r/2 with a being 
the Pauli matrices; g is the g-factor of the electron, which is chosen to be g = 2; /ib is 
the Bohr magneton of the electron; and Ai is the hyperhne coupling tensor, a 3x3 matrix. 
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As proposed in Ref. 05. we model the radical-pair dynamics with a Liouville equa¬ 
tion |162j . 


pit) = ~ H,p{t)\ 

- y _ y {Q T ’p(t)} • ( 23 ) 

In Eq. (125]) . H is the Hamiltonian of the system; Q s is the singlet projection operator, i.e. 
Q s = |5')(5'|, and Q 1 = |T + )(T + | + |T 0 )(T 0 |+ | T_)(T_\ is the triplet projection operator, 
where \S) is the singlet state and (|T + ), |T 0 ), |T_)) are the triplet states [ 163] ; p(t) is the 
density matrix for the system; ks and hr are the decay rates for the singlet state and 
triplet states, respectively. 

Under the basic scheme, we assume that the initial state of the radical pair is a per¬ 
fect singlet state, |S) = ^ 75 (IT4') — |4T))- Therefore, the initial condition for the density 
matrix is; p( 0 ) = \In <S> Q s , where the electron spins are in their singlet states, and 
nuclear spins are in a completely mixed state, which is a 4x4 identity matrix. Assuming 
that the rate is independent of spin, the decay rates for the singlet and triplet should be 
the same [41], kg = hr = k = l/zs -1 , i.e., k is the recombination rate for both the singlet 
and triplet states. The external weak magnetic held, B , representing the Earth’s magnetic 
held in Eq. (]22]h depends on the angles, 6 and p, with respect to the reference frame of 
the immobilized radical pair, i.e., B = R 0 (sin 6 cos 99 , sin 9 sin ip, cos 0), where B 0 = 0.5G 
is the magnitude of the local geomagnetic held. Without losing the essential physics, ip 
can be assumed to be 0. 


Since the radical pair must be very sensitive to different alignments of the magnetic held, 







Figure 4. Entanglements for different angles. Results are produced by assuming the decay 
rate to be 1/is -1 . All curves are almost identical. In the geomagnetic field, entanglement does 
not change with orientation. 


it is necessary to assume that the hyperhne coupling tensors in Eq. (|22|) are anisotropic. 
However, for the sake of simplicity, we employ the hyperhne coupling as anisotropic for 
one radical, and as isotropic for the other [U], i.e., 

n n \ 


( 10G 0 0 ^ 


A\ — 


\ 


0 10G 0 
0 0 0 


) ^2 — 


V 


5G 0 0 

0 5G 0 
0 0 5G 


Using the parameters defined above, we calculate one of the properties of the avian com¬ 
pass, depending only on the inclination but not on the polarity. As one can see in Fig. El 
the triplet yields are symmetric around 90° since the hyperhne coupling tensors are sym¬ 
metric. Consequently, the radical pair mechanism cannot distinguish between magnetic 
helds that are oppositely directed but have the same magnitude. 

Also, we investigated whether the entanglement, measured by negativity, is angle-dependent. 
While, using the suggested hyperhne coupling tensor in Ref. [41], the calculation gives 
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us the surprising result shown in Fig. [4j Namely, the dynamics of entanglement does 
not change with angle, i.e., the entanglement is not sensitive to the angle between the 
z-axis of the radical pair and the Earth’s magnetic held under these parameters. The 
main reason may be that the fast exponential decay rate hides the effect of the hyperhne 
coupling interaction. Therefore, when enlarging the hyperhne coupling, the difference of 
entanglement between angles will be observable. Even though the entanglement of the 
radical pair cannot directly provide directional information, this does not mean that en¬ 
tanglement is not involved in avian navigation. There might be indirect mechanisms for 
birds to utilize entanglement. 


The results of our calculations shown in Fig. [4l demonstrate that the dynamics of entan¬ 
glement is almost the same for all angles when symmetric hyperhne tensors are involved. 
This raises the following question: “What will happen if there is an asymmetric hyperhne 
tensor?” Although hyperhne tensors of organic radicals are usually symmetric, we can 
examine a few asymmetrical cases to try to find the possible asymmetry effects of the 
hyperhne coupling. The radicals with asymmetric hyperhne tensors are likely to have 
fast spin relaxation and not suitable for the purpose of a geomagnetic sensor. Therefore, 
our study here is only out of theoretical curiosity to find out what will happen to the 
entanglement when the hyperhne tensors are not symmetric. The asymmetric hyperhne 
tensors we examine are: 


A\ = 


10G 0 0 

0 10G 0 

0 0 4G 


\ / 

,#2 = 


V 


5G 5G 0 
0 5G 0 
0 0 5G 


\ 
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Figure 5. Entanglements for four angles under the hyperfine coupling tensors, A b i and A?. 


and 




^ 0 0 0 ^ 
0 0 0 
v 0 0 AG J 



( 0 5G 0 ^ 
0 0 0 
v 0 0 0 ) 


From Fig. 0 we can easily see that the hyperfine coupling tensor pair of gives an 
intriguing result. Namely, the dynamics of the entanglement is clearly dependent on the 
system’s orientation. 


3.2 Two-Stage Scheme 

For further study, we modify our model based on Ref. |98j. The radical pair reaction 
scheme has two stages (see Fig. (HD- The initial radical pair [FAD*“TrpH* + ] is formed by 
light-induced electron transfer, followed by the protonation and deprotonation, forming 
a secondary radical pair [FADH'Trp*]. This two-stage scheme is shown in Fig. [dj Both 
radical pairs are affected not only by the external magnetic field, but also by their sur¬ 
rounding nuclei. Respectively, the Hamiltonians of the initial and secondary radical pair 
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are, 


2 


Hi — g/JB Si ■ 

i= 1 

\B + An ■ Iij , 

(24) 

2 

H 2 = gpB Si ■ 

2—1 

( B + A2i ■ Ii'j . 

(25) 
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[FAD'TrpH‘ + ]^-[FAD‘TrpH ' + ] 3 
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f 


FAD+TrpH 


Triplet Products 


Figure 6. The reaction scheme of the radical pair mechanism in cryptochrome, kb = r b 1 and 
kf = rj 1 are the first-order rate constants for recombination of the initial radical pair and 
formation of the secondary pair from the initial one, respectively. k s = is rate constant for 
the decay of the secondary pair. The green two-headed arrows indicate the interconversion of 
the singlet and triplet states of the radical pairs. 


In Eq. (EMI) and Eq. (j25|) . Si is the unpaired electron spin of the radical pairs, and J* 
is the nuclear spin of nitrogen in the pairs. We calculate the hyperfine coupling tensors 
(Table. 1), A,-, using Gaussian09 with UB3LYP/EPR-II. For simplicity, in our subsequent 
calculations, we only use one of the hyperfine coupling tensors within each molecule, since 
additional nuclear spins have little effect on the yield curves E3- Also, because the 
electron is located near the nitrogen atoms, and the couplings between the electron and 
the nitrogen atoms are stronger than the couplings to other near-by hydrogen atoms, we 
choose, for our subsequent calculations, the hyperfine coupling tensors associated with the 
nitrogen atoms in each molecule. B is the weak external geomagnetic field. B depends on 
the angles, 6 and p, with respect to the reference frame of the immobilized radical pair, 
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i.e., B = Bo(sin6 cosip, sin 0 sin (p, cos6). We can choose the x-axis so that the azimuthal 
angle, ip, is 0. The constants, g and /is, are the //-factor and the Bohr magneton of the 
electron, respectively. 


Table. 1: The hyperfine coupling tensors of some atoms. 


Moleclule 

Atom 

Isotropic (G) 

Anisotropic (G) 

Principal Axis 

(Angstrom) 

TrpH ,+ 

N 

12.864 

-7.154 

0.63 0.70 -0.34 




-7.051 

0.73 0.63 0.04 




14.205 

-0.26 0.22 0.94 


H 

-3.054 

-1.478 

0.50 0.86 -0.12 




-0.977 

-0.14 0.20 0.96 




2.454 

0.85 -0.47 0.22 

FAD*” 

N 

2.339 

-5.392 

0.61 0.79 0.00 




-5.353 

0.79 0.61 0.00 




10.745 

0.00 0.00 1.00 

Trp* 

N 

8.393 

-9.708 

-0.23 0.97 0.10 




-9.539 

0.96 0.24 -0.14 




19.247 

0.15 -0.07 0.99 


H 

-5.888 

-2.458 

0.62 0.78 -0.09 




-0.792 

0.21 -0.06 0.98 




1.320 

0.75 -0.63 -0.20 

FADH* 

N 

2.015 

-4.815 

0.79 0.61 0.00 




-4.702 

0.61 0.79 0.00 




9.517 

0.00 0.00 1.00 


H 

-3.054 

-1.478 

0.50 0.86 -0.12 




-0.977 

-0.14 0.20 0.96 




2.454 

0.85 -0.47 0.22 
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The time evolution of the corresponding spin system is described through a modified 
stochastic Liouville equation [4T1I1641 - I168j . For this purpose, we denote the density matrix 
corresponding to the states of the radical pair [FAD*”TrpH* + ] as p\ and the density matrix 
corresponding to the states of the radical pair [FADH'Trp*] as [)•>. Each density matrix 
follows a stochastic Liouville equation that describes the spin motion and also takes into 
account the transition into and out of a particular state from or into other states, as 
illustrated in FigjG] Therefore, the dynamics of the radical pairs in the two-stage scheme 
is governed by the following coupled Liouville equations: 


dpi(t) 

dt 


y {< 2 5 ,PiW} - y {Q T ,Pi(t)} 
y {Q S ,Pi(t)} , 


(26) 


dpz(t) 

dt 


- ^[H 2 ,p 2 (t)] 

+ y {Q S ,Pi(t)} + y {Q T ,Pi(t)} 
~ y {Q S ,P2(t)} ~ y {Q T :P2(t)}, 


(27) 


where H i and H 2 are the Hamiltonians of the two radical pairs given in Eqs. (12T)i 
and (l25]k Q s , as defined before, is the singlet projection operator, Q s = |S')(S'|, and 
Q T = |T + )(T + | + |To)(To| + |XL)(XL| is the triplet projection operator, where (S') is the 
singlet state, and (|X + ), |X 0 ), |X_)) are the triplet states; and all of the decay rates are 
indicated in Fig. [6j In addition, the initial state of the pair [FAD*”TrpH ,+ ] is assumed to 
be in the singlet state, |S) = -^(ITI) - IIT))> while the pair [FADH'Trp*] is not produced 
initially. In other words, Pi(0) = 1 1^ <S> Q s , where the electron spins are in the singlet 
states, and nuclear spins are in thermal equilibrium, a completely mixed state, which is 
a 9x9 identity matrix, and ^ 2 ( 0 ) = 0. 
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We consider the product formed by the radical pair [FADH'Trp*] in the triplet state 
as the signal product, whose yield is defined as: = k s J 0 °° Tr[Q T p{t)]dt [ 881I16511T69] . 

where Q T = \T)(T\, and |T) = |T+> + \T 0 ) + \T_). 



Figure 7. Entanglement of the initial radical pair [FAD*~TrpH* + ] as a function of the 
external fields for four polar angles, 6 = 0°(black), 30°(red), 60°(blue), 90°(green). Since the 
entanglement of the initial pair [FAD*~TrpH* + ] is compressed within 0.1 /rs, the timescale of 
the graph is from 0 to 0.1 /us. The other graphs range from 0 to 0.8 ^s. And the entanglements 
of the initial pair at 0°,30°,60° differ after 0.3 /rs. 


The entanglement is believed to play an important role in many systems | I1 54111 7411176] . 
including the chemical compass in birds. As mentioned before, we use negativity as the 
metric of entanglement. However, for the two-stage scheme, the secondary radical pair 
[FADH'Trp*] barely has any entanglement between the two unpaired electrons, since the 
chemical reaction (protonation and deprotonation) has destroyed the entanglement be¬ 
tween them in the preceding radical pair [FAD*~TrpH* + ], The unpaired electrons in the 
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initial radical pair show a robust entanglement. Fig. [7] shows the entanglement of the 
initial radical pair [FAD* _ TrpH* + ] for four polar angles, 6. Also, the dynamics of the 
entanglement is clearly dependent on the angles, which is very different from the results 
in the one-stage case [lfTO] . However, the entanglements at the angles, 0°, 30° and 60°, 
are nearly the same for the first 0.1/rs, while the entanglement at 90° is very different 
from others. At 90°, the entanglement lasts for 0.1/zs, which is long enough for electrons 
to transfer between different molecules mg. 


The triplet yields of FADH+Trp The Triple Yields of FAD + TrpH 
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Figure 8. The triplet yields of radical pairs [FADH'Trp*] and [FAD*~TrpH* + ] as a 
function of angles. The figures show that the yield of the pair [FAD* _ TrpH* + ] is almost 
zero. That means, after a relative long time, the radical pair [FAD* _ TrpH* + ] has 
converted to the pair [FADH'Trp*] via chemical reactions. Thus, the triplet yield of 
[FAD*”TrpH* + ] is almost zero. 


Fig. [8]shows the yields of [FAD*”TrpH* + ] and [FADH'Trp*]. We can tell that after a rel¬ 
atively long time, the triplet yield of [FAD*~TrpH* + ] is almost zero, which demonstrates 
that the pair of [FAD* TrpH* + ] has transferred to [FADH'Trp*] via some chemical reac¬ 
tions. Basically, following the scheme we used, these results of yields can be important, 
and the yield of [FADH'Trp*] can be seen as the signal for birds. Around 90° (80°-90°), 
the derivative of yields with respect to the angle seems to be larger than for the other 
angles, which indicates that the birds are more sensitive when they are heading north. 
This could be a good sign, because it may give birds the cue of direction. 
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0 (degrees) 


Figure 9. The magnetic sensitivity of the chemical compass as a function of the angle. The 
sensitivity is defined as d&x/dB, in T _1 . There is a rapid increase in this sensitivity between 
80° and 90°. The sensitivity under geomagnetic field is of the order of 1CF 3 , requiring a strong 
magnification mechanism in birds to utilize it. 
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We now focus on the magnetic sensitivity of the avian compass, which is defined as 
chhr /dB (T" 1 ) |170UT71 ]. Fig. |9] shows that the sensitivities around 0° and 90° are sim¬ 
ilar and also larger than for most other angles, which could indicate that the birds can 
detect the directions of meridians and parallels if they use the intensity of the magnetic 
held for navigation, since the yield-based compass is most sensitive along these two di¬ 
rections. Another property that attracted our attention is that the sensitivity’s slope is 
significantly larger between 80° and 90° than that of the other sections of the curve. This 
property of increased sensitivity may imply that it is easier for birds to detect the direc¬ 
tion of magnetic parallels than that of magnetic meridians. Since the yield-based compass 
is very sensitive to the change of intensities, we can also expect that it is easier for birds 
to detect the change of the held intensities when the polar angle is around 90°. This 
capability can enable birds to migrate along the direction of the gradient of intensities of 
the geomagnetic helds. 

Furthermore, we explore the magnetic sensitivity as a function of the intensities of the 
magnetic held for several polar angles 9, in Fig. |T0l Fig. flOl shows two types of curves. 
The hrst pattern is observed for 85° and 90°, in which the sensitivities monotonically de¬ 
crease as the external helds increase. In this situation, the sensitivities are much higher for 
very weak magnetic helds, less than 0.25G, than those in the normal range of the geomag¬ 
netic helds, from 0.25G to 0.65G cm The sensitivities fall into the normal range in the 
geomagnetic helds, similar to other angles. The other pattern occurs for 0°, 30° and 60°, 
and the sensitivities increase initially, and then decrease as the external helds increase. 
In this situation, the maxima of the curves move rightwards and downwards as the polar 
angles increase. Combining these two situations (Fig. fTOlh we observe the properties of 
the chemical compass mentioned before, namely that compass is most magnetically sensi¬ 
tive around 0° and 90° at the geomagnetic helds. However, above 0.35G, all sensitivities 
decrease as the helds’ intensities increase. This may explain why some species of birds 
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Figure 10. (Color online) Intensity dependence of the magnetic sensitivity. This graph 
shows the sensitivity as a function of external magnetic field, B, as a function of polar angle, 9, 
between the z-axis of the radical pair and the magnetic field, i.e., 9 = 0 °(black), 30 °(red), 

60°(blue), 85°(pink), 90 °(green). The interior graph magnifies the data for angles 
6 = 0° (black), 30° (red), 60° (blue). 
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lose their ability to orient themselves in higher magnetic fields [HI 175] . Also, since the 
sensitivity is not zero, after extended exposure to unnatural magnetic fields the birds may 
adapt to the decreased sensitivity, so that they are able to regain the ability to orient [173] . 

All of the above results can provide us with a basic picture of the radical pair mech¬ 
anism. 


4 Conclusions 

The radical pair mechanism is a promising hypothesis to explain the mystery of the navi¬ 
gation of birds. This theoretical study has demonstrated the role of weak magnetic fields 
play in the product yields of the radical pairs. In addition, this type of study has inspired 
scientists to design highly effective devices to detect weak magnetic fields and to use the 
geomagnetic fields to navigate. 

The anisotropic hyperfine coupling between the electron spins and the surrounding nu¬ 
clear spins can play a crucial role in avian magnetoreception. The hyperfine coupling can 
affect not only the product yields but also the entanglement of the electron spin states. 
By involving more nuclear spins one can greatly enhance the quantum entanglement [ 478j . 
Additionally, mimicking this anisotropic magnetic environment can be very useful for cre¬ 
ating detectors of weak magnetic fields. 

By studying the role of intensity of the magnetic field in avian navigation, we find that 
birds could be able to detect the change of the intensity of geomagnetic fields and the 
approximate direction of parallels instead of sensing the exact direction. However, the 
mechanism in which birds can utilize the signal remains unknown at this time. 


40 





Acknowledgements 

S.K. and Y.Z. would like to thank the NSF Center for Quantum Information for Quantum 
Chemistry (QIQC), Award No. CHE-1037992, for financial support. The work by G.P.B. 
was carried out under the auspices of the National Nuclear Security Administration of 
the U.S. Department of Energy at Los Alamos National Laboratory under Contract No. 
DE-AC52-06NA25396. 

References 

[1] H. Hayashi, Introduction to Dynamic Spin Chemistry. Magnetic Field Effects on 
Chemical and Biochemical Reactions. World Scientific, 2004 

[2] S. Nagakura, H. Hayashi and T. Azumi, Dynamic Spin Chemistry. Magnetic Con¬ 
trols and Spin Dynamics of Chemical Reactions. A KODANSHA LTD. and WILEY 
& SONS, INC., 1998 

[3] R. Z. Sagdeev, Yu. N. Molin, K. M. Salikhov, T. V. Leshina, M. A. Kamha and S. 
M. Shein: Effects of magnetic held on chemical reactions. Org. Magn. Res., 5, 603 
(1973) 

[4] B. Brocklehurst, R. S. Dixon, E. M. Gardy, V. J. Ropata, M. J. Quinn, A. Shin and 
F. P. Sargent: The effect of a magnetic held on the singlet/triplet ratio in geminate 
ion recombination. Chem. Phys. Lett., 28, 361 (1974) 

[5] K. Schulten, H. Staerk, A. Weller, H. J. Werner and B. Nickel: Magnetic held 
dependence of the geminate recombination of radical ion pairs in polar solvents. Z. 
Phys. Chem. (Frankfurt/Main), 101, 371 (1976) 


41 


[ 6 ] Y. Tanimoto, H. Hayashi, S. Nagakura, H. Sakuragi and K. Tokumaru: The external 
magnetic field effect on the singlet sensitized photolysis of dibenzoyl peroxide. Chem. 
Phys. Lett., 41, 267 (1976) 

[7] N. Hata: The effect of external magnetic field on the photochemical reaction of 
isoquinoline N-oxide. Chem. Lett., 547 (1976) 

[ 8 ] Z. Schulten and K. Schulten: The generation, diffusion, spin motion, and recombi¬ 
nation of radical pairs in solution in the nanosecond time domain. J. Chem. Phys., 
66 , 4616 (1977) 

[9] A. M. Lewis, D. E. Manolopoulos and P. J. Hore: Asymmetric recombination and 
electron spin relaxation in the semiclassical theory of radical pair reactions. J. Chem. 
Phys., 141, 044111 (2014) 

[10] H. Hayashi and S. Nagakura: The Theoretical Study of External Magnetic Field 
Effect on Chemical Reactions in Solution. Bull. Chem. Soc. Jpn., 51, 2862 (1978) 

[11] R. Kaptein: Chemically induced dynamic nuclear polarization. VIII. Spin dynamics 
and diffusion of radical pairs. J. Am. Chem. Soc., 94, 6251 (1972) 

[12] R. Kaptein and J. A. Den Hollander: Chemically induced dynamic nuclear polar¬ 
ization. X. Magnetic field dependence. J. Am. Chem. Soc., 94, 6269 (1972) 

[13] A. Weller, F. Nolting and H. Staerk: A quantitative interpretation of the magnetic 
field effect on hyperfine-coupling-induced triplet fromation from radical ion pairs. 
Chem. Phys. Lett., 96, 24 (1983) 

[14] E. Schodinger, What is life? Cambridge Univ. Press, 1992 

[15] M. Mohseni, Y. Omar, G. Engel and M. Plcnio, Quantum effects in biology. Cam¬ 
bridge University Press (2014). 


42 


[16] N. Lambert, Y. Chen, Y. Cheng, C. Li, G. Chen and F. Nori: Quantum biology. 
Nature Physics, 9, 10 (2013) 

[17] P. Ball: Physics of life: The dawn of quantum biology. Nature (London), 474, 272 

( 2011 ) 

[18] H. C. Longuet-Higgins: Quantum Mechanics and Biology. Biophy. J., 2, 207 (1962) 

[19] A. Shabani, M. Mohseni, S. Jang, A. Ishizaki, M. B. Plcnio, P. Rebentrost, A. 
Aspuru-Guzik, J. Cao, S. Lloyd and R. Silbey: Open quantum system approaches 
to biological systems. In Quantum effects in biology ed. M. Mohseni, Y. Omar, G. 
Engel and M. Plcnio, Cambridge University Press (2014). 

[20] G. Panitchayangkoon, D. Hayes, K. A. Fransted, J. R. Caram, E. Harel, J. Wen, R. 
E. Blankenship and G. S. Engel: Long-lived quantum coherence in photosynthetic 
complexes at physiological temperature. Proc. Natl Acad. Sci., 107, 12766 (2010) 

[21] P. Rebentrost, M. Mohseni, I. Kassal, S. Lloyd and A. Aspuru-Guzik: Environment- 
assisted quantum transport. New J. Phys., 11, 033003 (2009) 

[22] M. Mohseni, P. Rebentrost, S. Lloyd and A. Aspuru-Guzik: Environment-assisted 
quantum walks in photosynthetic energy transfer. J. Chem. Phys., 129, 174106 
(2008) 

[23] A. Ishizaki and G. R. Fleming: Theoretical examination of quantum coherence in 
a photosynthetic system at physiological temperature. Proc. Natl Acad. Sci., 106, 
17255 (2009) 

[24] C. Li, N. Lambert, Y. Chen, G. Chen and F. Nori: Witnessing Quantum Coherence: 
from solid-state to biological systems. Sci. Rep., 2, 885 ( 2012 ) 


43 


[25] A. S. Fokas, D. J. Cole and A. W. Chin: Constrained geometric dynamics of the 
Fenna-Matthews-Olson complex: the role of correlated motion in reducing uncer¬ 
tainty in excitation energy transfer. Photosynth. Res., 122, 275 (2014) 

[26] F. Caruso, A.W. Chin, A. Datta, S.F. Huelga and M.B. Plenio: Highly efficient 
energy excitation transfer in light-harvesting complexes: The fundamental role of 
noise-assisted transport. J. Chem. Phys., 131, 105106 (2009) 

[27] G. S. Engel, T. R. Calhoun, E. L. Read, T. K. Aim, T. Mancal, Y-C. Cheng, R. 
E. Blankenship and G. R. Fleming: Evidence for wavelike energy transfer through 
quantum coherence in photosynthetic complexes. Nature, 446, 782 (2007) 

[28] D. Hayes and G. S. Engel: Extracting the excitonic Hamiltonian of the Fenna- 
Matthews-Olson complex using three-dimensional electronic spectroscopy. Biophys. 
J., 100, 2043 (2011) 

[29] G. S. Engel: Quantum coherence in photosynthesis. Chemistry Procedia, 3, 222 

( 2011 ) 

[30] G. Panitchayangkoon, D. V. Voronine, D. Abramavicius, J. R. Caram, N. H. C. 
Lewis, S. Mukamel and G.S. Engel: Direct evidence of quantum transport in pho¬ 
tosynthetic light-harvesting complexes. Proc. Natl Acad. Sci., 108, 20908 (2011) 

[31] J.R. Caram and G.S. Engel: Extracting dynamics of excitonic coherences in con¬ 
gested spectra of photosynthetic light harvesting antenna complexes. Farad. Dis¬ 
cuss., 153, 93 (2011) 

[32] A. F. Fidler, J. R. Caram, D. Hayes and G.S. Engel: Toward a coherent picture of 
excitonic coherence in the Fenna?Matthews?01son complex. J. Phys. B, 45, 154013 
( 2012 ) 


44 


[33] K. M. Pelzer, T. Can, S. K. Gray, D. K. Morr and G. S. Engel: Coherent transport 
and energy flow patterns in photosynthesis under incoherent excitation. J. Phys. 
Chem. B, 118, 2693 (2014) 

[34] M. Sarovar, A. Ishizaki, G. R. Fleming and K. B. Whaley: Quantum entanglement 
in photosynthetic light-harvesting complexes. Nature Phys., 6, 462 (2010) 

[35] S. Hoyer, M. Sarovar and K. B. Whaley: Limits of quantum speedup in photosyn¬ 
thetic light harvesting. New J. Phys., 12, 065041 (2010) 

[36] M. Sarovar, Y.-C. Cheng and K. B. Whaley: Environmental correlation effects on 
excitation energy transfer in photosynthetic light harvesting. Phys. Rev. E, 83, 
011906 (2011) 

[37] S. Hoyer, F. Caruso, S. Montangero, M. Sarovar, T. Calarco, M. B. Plenio and K. 
B. Whaley: Realistic and verifiable coherent control of excitonic states in a light¬ 
harvesting complex. New J. Phys., 16, 045007 (2014) 

[38] J. Zhu, S. Kais, A. Aspuru-Guzik, S. Rodriques, B. Brock and P. J. Love: Multipar¬ 
tite quantum entanglement evolution in photosynthetic complexes. J. Chem. Phys., 
137, 074112 (2012) 

[39] F. Caruso, A. W. Chin, A. Datta, S. F. Huelga and M. B. Plenio: Entanglement 
and entangling power of the dynamics in light-harvesting complexes. Phys. Rev. A, 
81, 062346 (2010) 

[40] F. Caycedo-Soler, A. W. Chin, J. Almeida, S. F. Huelga and M. B. Plenio: The 
nature of the low energy band of the Fenna-Matthews-Olson complex: vibronic 
signatures. J. Chem. Phys., 136, 155102 (2012) 


45 


[41] T. Ritz, S. Adem and K. Schulten: A model for photoreceptor-based magnetore¬ 
ception in birds. Biophy. J., 78, 707 (2000) 

[42] W. Wiltschko and R. Wiltschko: Magnetic Compass of European Robins. Science, 
176, 62 (1972) 

[43] K. Schulten, C. E. Swenberg and A. Weller: A biomagnetic sensory mechanism 
based on magnetic field modulated coherent electron spin motion. Z. Phys. Chem., 
Ill, 1 (1978) 

[44] K. Schulten: Magnetic field effects in chemistry and biology. Festkorperprobleme, 
22, 61 (1982) 

[45] U. E. Steiner and T. Ulrich: Magnetic field effects in chemical kinetics and related 
phenomena. Chem. Rev., 89, 51 (1989) 

[46] W. Wiltschko, K. Stapput, P. Thalau and R. Wiltschko: Avian magnetic com¬ 
pass: fast adjustment to intensities outside the normal functional window. Natur- 
wissenschaften, 93, 300 (2006) 

[47] C. T. Rodgers and P. J. Hore: Chemical magnetoreception in birds: the radical pair 
mechanism. Proc. Natl Acad. Sci., 106, 353 (2009) 

[48] T. Ritz: Quantum effects in biology: Bird navigation. Proc. Chem., 3, 262 (2011) 

[49] P. Thalau, T. Ritz, H. Burda, R. E. Wegner and R. Wiltschko: The magnetic 
compass mechanisms of birds and rodents are based on different physical principles. 
J. R. Soc. Interface, 3, 583 (2006) 

[50] W. Wiltschko and R. Wiltschko: Magnetoreception in birds: two receptors for two 
different tasks. J. Ornithol, 148, S61 (2007) 


46 


[51] N. Lambert, S. De Liberato, C. Emary and F. Nori: Radical-pair model of magneto 
reception with spin-orbit coupling. New J. Phys., 15, 083024 (2013) 

[52] S. F. Huelga and M. B. Plenio: Vibriations, Quanta and Biology. Contemp. Phys., 
54, 181 (2013) 

[53] S. R. T. Neil, J. Li, D. M.W. Sheppard, J. Storey, K. Maeda, K. B. Henbest, P. 
J. Hore, C. R. Timmel and S. R. Mackenzie: Broadband cavity-enhanced detection 
of magnetic field effects in chemical models of a cryptochrome magnetoreceptor. J. 
Phys. Chem. B ., 118, 4177 (2014) 

[54] I. A. Solov’yov, T. Ritz, K. Schulten and P. J. Hore: A chemical compass for bird 
navigation. In Quantum effects in biology ed. M. Mohseni, Y. Omar, G. Engel and 
M. Plenio, Cambridge University Press (2014). 

[55] F. Cintolesi, T. Ritz, C.W.M. Kay, C. R. Timmel and P. J. Hore: Anisotropic 
recombination of an immobilized photoinduced radical pair in a 50-/iT magnetic 
field: a model avian photomagnetoreceptor. Chem. Phys., 294, 385 (2003) 

[56] L. Turin: A Spectroscopic Mechanism for Primary Olfactory Reception. Chem. 
Senses, 21, 773 (1996) 

[57] J. C. Brookes, F. Hartoutsiou, A. P. Horsfield and A. M. Stoneham: Could humans 
recognize odor by phonon assisted tunneling? Phys. Rev. Lett., 98, 038101 (2007) 

[58] V. P. Bingman, K. P. Able and P. Kerlinger: Wind drift, compensation, and the 
use of landmarks by nocturnal bird migrants. Anim. Behav., 30, 49 (1982) 

[59] T. C. Williams, J. M. Williams, P. G. Williams and P. Stokstad: Bird migration 
through a mountain pass studied with high resolution radar, ceilometers, and census. 
The Auk, 118, 389 (2001) 


47 


[60] S. T. Emlen: Migratory orientation in the Indigo Bunting Passerina cyanea. Part 
I: evidence for use of celestial cues. The Auk, 84, 309 (1967) 

[61] S. T. Emlen: Migratory orientation in the Indigo Bunting Passerina cyanea. Part 
II: Mechanism of celestial orientation. The Auk, 84, 463 (1967) 

[62] S. T. Emlen: The development of migratory orientation in young Indigo Buntings. 
Living Bird , 8, 113 (1969) 

[63] S. T. Emlen: Celestial rotation: its importance in the development of migratory 
orientation. Science, 170, 1198 (1970) 

[64] P. Weindler, B. Mathias and W. Wiltschko: The direction of celestial rotation 
influences the development of stellar orientation in young Garden Warblers ( Sylvia 
Borin). J. Exp. Biol., 200, 2017 (1997) 

[65] R. Wiltschko, C. Haugh, M. Walker and W. Wiltschko: Pigeon homing: sun compass 
use in the southern hemisphere. Behav. Ecol. Sociobiol. , 43, 297 (1998) 

[66] S. J.Duff, L. A. Brownlie, D. F. Sherry, and M. Sangster: Sun compass and landmark 
orientation by Black-capped Chickadees ( Parus atricapillus). J. Exp. Psychol. Anim. 
Behav. Process, 24, 243 (1998) 

[67] R. Wiltschko, M. Walker, and W. Wiltschko: Sun-compass orientation in Homing 
Pigeons: compensation for different rates of change in azimuth? J. Exp. Biol, 203 
889 (2000) 

[68] T. Ritz, P. Tlialau, J. B. Phillips, R. Wiltschko and W. Wiltschko: Resonance 
effects indicate a radical-pair mechanism for avian magnetic compass. Nature, 429, 
177 (2004) 


48 


[69] M. Roy, S. K. Meena, S. K. Singh, N. K. Sethy, K. Bhargava, S. Sarkar and M. 
Das: Presence of stable carbon centric free radicals and ferromagnetic elements in 
the antennae and the wings of nocturnal silk moth: A magnetic nanostructure for 
magneto sensing. Mater. Express, 3, 43 (2013) 

[70] U. Munro, J. A. Munro and J. B. Phillips: Evidence for a magnetite-based naviga¬ 
tional ’’map” in birds. Naturwissenschaften, 84, 26 (1997) 

[71] J. L. Kirschvink, M. M. Walker and C. E. Diebel: Magnetite-based magnetorecep¬ 
tion. Curr. Opin. Neurobiol., 11, 462 (2001) 

[72] C. D. Treiber et al.: Clusters of iron-rich cells in the upper beak of pigeons are 
macrophages not magnetosensitive neurons. Nature, 484, 367 (2012) 

[73] W. Wiltschko, R. Freire, U. Munro, T. Ritz, L. Rogers, P. Tlialau and R. Wiltschko: 
The magnetic compass of domestic chickens, Gallus gallus. J. Exp. Biology, 210, 
2300 (2007) 

[74] W. Wiltschko and R. Wiltschko: Magnetic orientation and magnetoreception in 
birds and other animals. J. Comp. Physiol. A, 191, 675 (2005) 

[75] W. Wiltschko and R. Wiltschko: Magnetic Compass of European Robins. Science, 
176, 62 (1972) 

[76] R. Wiltschko and W. Wiltschko: Avian magnetic compass: its functional properties 
and physical basis. Curr. Zool., 56, 265 (2010) 

[77] A. F. Davila, G. Fleissner, M. Winklhofer and N. Petersen: A new model fro a mag¬ 
neto receptor in homing pigeons based on interacting clusters of superparamagnetic 
magnetite. Phys. Chem. Earth, 28, 647 (2003) 


49 


[78] W. Wiltschko, U. Munro, H. Ford and R. Wiltschko: Red light disrupts magnetic 
orientation of migratory birds. Nature, 364, 525 (1993) 

[79] W. Wiltschko and R. Wiltschko: Migratory orientation of European Robins is af¬ 
fected by the wavelength of light as well as by a magnetic pulse. J. Comp. Physiol. 
A, 177, 363 (1995) 

[80] W. Wiltschko and R. Wiltschko: Pigeon Homing: Effect of Various Wavelengths of 
Light During Displacement. Naturwissenschaften, 85, 164 (1998) 

[81] W. Wiltschko and R. Wiltschko: The effect of yellow and blue light on magnetic 
compass orientation in European robins, Erithacus rubecula. J. Comp. Physiol. A, 
184, 295 (1999) 

[82] R. Wiltschko, T. Ritz, K. Stapput, P. Thalau and W. Wiltschko: Two different 
types of light-dependent responses to magnetic fields in birds, curr. biol, 15, 1518 
(2005) 

[83] K. Stapput, P. Thalau, R. Wiltschko and W. Wiltschko: Orientation of birds in 
total darkness. Curr. Biol., 18, 602 (2008) 

[84] R. Wiltschko, S. Denzau, D. Gehring, P. Thalau, and W. Wiltschko: Magnetic 
orientation of migratory robins, Erithacus rubecula, under long-wavelength light. J. 
Exp. Biology, 214, 3096 (2011) 

[85] T. Ritz, R. Wiltschko, P. J. Hore, C. T. Rodgers, K. Stapput, P. Thalau, C. R. 
Timmcl and W. Wiltschko: Magnetic compass of birds is based on a molecule with 
optimal directional sensitivity. Biophys. J., 96, 3451 (2009) 


50 


[86] J. C. S. Lau, N. Wagner-Rundell, C. T. Rodgers, N. J. B. Green and P. J. Hore: Ef¬ 
fects of disorder and motion in a radical pair magneroreceptor. J. R. Soc. Interface, 
7 , S257 (2010) 

[87] C. A. Dodson, P. J. Hore and M. I. Wallace: A radical sense of direction: signalling 
and mechanism in cryptochrome magneroreception. Trends Biochem. Sci., 38 , 435 
(2013) 

[88] C. R. Tinnncl, U. Till, B. Brocklehurst, K. A . McLauchlan and P. J. Hore: Effects 
of weak magnetic fields on free radical recombination reactions. Mol. Phys., 95 , 71 
(1998). 

[89] W. Wiltschko, J. Traudt, O. Giintiirkiin, H. Prior and R. Wiltschko: Lateralisation 
of magnetic compass orientation in a migratory bird. Nature, 419 , 467 (2002) 

[90] K. Stapput, O. Giintiirkiin, K. P. Hoffmann, R. Wiltschko and W. Wiltschko: Mag¬ 
netoreception of directional information in birds requires nondegraded vision. Curr. 
Biol, 20 , 1259 (2010) 

[91] H. Mouritsen, U. Janssen-Bienhold, M. Liedvogel, G. Feenders, J. Stalleicken, P. 
Dirks and R. Weiler: Cryptochromes and neuronal-activity markers colocalize in 
the retina of migratory birds during magnetic orientation. Proc. Natl Acad. Sci., 
101 , 14294 (2004) 

[92] T. Ritz, T. Yosliii, C. Helfrich-Foerster and M. Ahmad: Cryptochrome: A pho¬ 
toreceptor with the properties of a magnetoreceptor? Commun Integr Biol., 3, 24 
( 2010 ) 

[93] C. Nie/3ner, S. Denzau, K. Stapput, M. Ahmad, L. Peichl, W. Wiltschko and R. 
Wiltschko: Magnetoreception: activated cryptochrome la concurs with magnetic 
orientation in birds. J. R. Soc. Interface, 10 , 20130638 (2013) 


51 


[94] S.-R. Harris, K. B. Henbest, K. Maeda, J. R. Panned, C. R. Timmel, P. J. Hore 
and H. Okamoto: Effect of magnetic fields on cryptochrome-dependent responses in 
Arabidopsis thaliana. J. R. Soc., Interface , 6 , 1193 (2009) 

[95] I. A. Solov’yov, H. Mouritsen and K. Schnlten: Acuity of a cryptochrome and 
vision-based magnetoreception system in birds. Biophys. J., 99 , 40 (2010) 

[96] J. C. S. Lau, C. T. Rodgers and P. J. Hore: Compass magnetoreception in birds 
arising from photo-induced radical pairs in rotationally disordered cryptochromes. 
J. R. Soc. Interface, 9 , 3329 (2012) 

[97] I. A. Solov’yov and K. Schnlten: Reaction Kinetics and Mechanism of Magnetic 
Field Effects in Cryptochrome. J. Phys. Chem. B, 116 , 1089 (2012) 

[98] K. Maeda, A. J. Robinson, K. B. Henbest, H. J. Hogben, T. Biskup, M. Ahmad, 
E. Schleicher, S. Weber, C. R. Timmel and P. J. Hore: Magnetically sensitive light- 
induced reactions in cryptochrome are consistent with its proposed role as a mag¬ 
netoreceptor. Proc. Natl Acad. Sci ., 109 , 4774 (2012) 

[99] K. Maeda, K. B. Henbest, F. Cintolesi, I. Kuprov, C. T. Rodgers, P. A. Liddell, D. 
Gust, C. R. Timmel and P. J. Hore: Chemical compass model of avian magnetore¬ 
ception. Nature, 453, 387 (2008) 

[100] A. Moller, S. Sagasser, W. Wiltschko and B. Schierwater: Retinal cryptochrome in 
a migratory passerine bird: a possible transducer for the avian magnetic compass. 
Naturwissenschaften, 91 , 585 (2004) 

[101] M. Liedvogel, K. Maeda, K. Henbest, E. Schleicher, T. Simon, C. R. Timmel, P. 
J. Hore and H. Mouritsen: Chemical Magnetoreception: Bird Cryptochrome la Is 
Excited by Blue Light and Forms Long-Lived Radical-Pairs. PLoS ONE, 2 , ell06 
(2007) 


52 


[102] T. Biskup, E. Schleicher, A. Okafuji, G. Link, K. Hitomi, E. D. Getzoff and S. 
Weber: Direct observation of a photoinduced radical pair in a cryptochrome blue- 
light photoreceptor. Angew. Chem., Int. Ed ., 48 , 404 (2009) 

[103] M. Liedvogel and H. Mouritsen: Cryptochromes-a potential magnetoreceptor: what 
do we know and what do we want to know? J. R. Soc. Interface, 7 , S147 (2010) 

[104] C. Nie/5ner, S. Denzau, J. C. Gross, L. Peichl, H.-J. Bischof, G. Fleissner, W. 
Wiltschko and R. Wiltschko: Avian Ultraviolet/Violet Cones Identified as Probable 
Magnetoreceptors. PLoS ONE, 6 , e20091 (2011) 

[105] M. Zapka, D. Heyers, C. M. Hein, S. Engels, N.-L. Schneider, J. Hans, S. Wcilcr, 
D. Dreyer, D. Kishkinev, J. M. Wild and H. Mouritsen: Visual but not trigeminal 
mediation of magnetic compass information in a migratory bird. Nature, 461 , 1274 
(2009) 

[106] H. Mouritsen, G. Feenders, M. Liedvogel, K. Wada and E. D. Jarvis: Night-vision 
brain area in migratory songbirds. Proc. Natl. Acad. Sci. U.S.A., 102 , 8339 (2005) 

[107] G. Feenders, M. Liedvogel, M. Rivas, M. Zapka, H. Horita, E. Hara, K. Wada, H. 
Mouritsen and E. D. Jarvis: Molecular Mapping of Movement-Associated Areas in 
the Avian Brain: A Motor Theory for Vocal Learning Origin. PLoS ONE, 3 , el768 
(2008) 

[108] D. Heyers, M. Manns, H. Luksch, O. Giintiirkiin and H. Mouritsen: A visual path¬ 
way links brain structures active during magnetic compass orientation in migratory 
birds. PLoS ONE, 2 , e937 (2007) 

[109] H. Mouritsen and P. J. Hore: The magnetic retina: light-dependent and trigeminal 
magneto reception in migratory birds. Curr. Opin. Neurobiol., 22 , 343 (2012) 


53 


[110] M. Ahmad, P. Galland, T. Ritz, R. Wiltschko and W. Wiltschko: Magnetic intensity 
affects crypochrome-dependent responses in Arabidopsis thaliana. Planta, 225 , 615 
(2007) 

[111] Y.-T. Kao, C. Tan, S.-H. Song, N. Oztiirk, J. Li, L. Wang, A. Sancar and D. Zhong: 
Ultrafast dynamics and anionic active states of the flavin cofactor in cryptochrome 
and photolyase. J. Am. Chem. Soc., 130 , 7695 (2008) 

[112] J.-P. Bouly, E. Schleicher, M. Dionisio-Sese, F. Vandenbussche, D. Van Der Straeten, 
N. Bakrirn, S. Meier, A. Batschauer, P. Galland, R. Bittl and M. Ahmad: Cryp¬ 
tochrome blue light photoreceptors are activated through interconversion of flavin 
redox states. J. Biol. Chem., 282 , 9383 (2007) 

[113] J. Shirdel, P. Zirak, A. Penzkofer, H. Breitkreuz and E. Wolf: Absorption and 
fluorescence spectroscopic characterisation of the circadian blue-light photoreceptor 
cryptochrome from Drosophila melanogaster (dCry). Chem. Phys., 352, 35 (2008) 

[114] N. ffoang, E. Schleicher, S. Kacprzak, J.P. Bouly, M. Picot, W. Wu, A. Berndt, E. 
Wolf, R.Bittl and M. Ahmad: Human and Dropsophila crypto chromes are light 
activated by flavin photo reduction in living cells. PLoS Biol., 6 , 1559 (2008) 

[115] K. B. Henbest, K. Maeda, P. J. Hore, M. Joshi, A. Bacher, R. Bittl, S. Weber, C. 
R. Timmel and E. Schleicher: Magnetic-field effect on the photoactivation reaction 
of Escherichia coli DNA photolyase. Proc. Natl Acad. Sci., 105 , 14395 (2008) 

[116] O. Efimova and P. J. Hore: Role of exchange and dipolar interaction in the radical 
pair model of the avian magnetic compass. Biophys. J., 94 , 1565 (2008) 

[117] I. A. Solov’yov and K. Schulten: Magnetoreception through cryptochrome may 
involve superoxide. Biophys. J., 96 , 4804 (2009) 


54 


[118] H. J. Hogben, O. Efimova, N. Wagner-Rundell, C. R. Timmel and P. J. Hore: 
Possible involvement of superoxide and dioxygen with cryptochrome in avian mag¬ 
neto reception: origin of Zeeman resonances observed by in vivo EPR spectroscopy. 
Chem. Phys. Lett., 480 , 118 (2009) 

[119] I. A. Solov’yov, D. Chandler, and K. Schulten: Magnetic Field Effects in Arabidopsis 
thaliana Cryptochrome-1. Biophys. J., 92 , 2711,(2007) 

[120] A. A. Lee, J. C. S. Lau, ff. J. Hogben, T. Biskup, D. R. Kattnig and P. J. Hore: 
Alternative radical pairs for cryptochrome-based magnetoreception. J. Roy. Soc. 
Interface., 11 , 20131063 (2014) 

[121] C. R. Timmel, F. Cintolesi, B. Brocklehurst and P. J. Hore: Model calculations of 
magnetic field effects on the recombination reactions of radicals with anisotropic 
hyperfine interactions. Chem. Phys. Lett., 334 , 387 (2001) 

[122] K. Maeda, C. J. Wedge, J. G. Storey, K. B. Henbest, P. A. Liddell, G. Kodis, D. 
Gust, P. J. Hore and C. R. Timmel: Spin-selective recombination kinetics of a model 
chemical magnetoreceptor. Chem. Commun., 47 , 6563 (2011) 

[123] K. Maeda, J. G. Storey, P. A. Liddell, D. Gust, P. J. Hore, C. J. Wedge and C. R. 
Timmel: Probing a chemical compass: novel variants of low-frequency reaction yield 
detected magnetic resonance. Phys. Chem. Chem. Phys., DOI: 10.1039/c4cp04095c 

[124] H. J. Hogben, T. Biskup and P. J. Hore: Entanglement and sources of magnetic 
anisotropy in radical pair-based avian magnetoreceptors. Phys. Rev. Lett., 109 , 
220501 (2012) 

[125] J. Cai, F. Caruso and M. B. Plcnio: Quantum limits for the magnetic sensitivity of 
a chemical compass. Phys. Rev. A, 85 , 040304 (2012) 


55 


[126] C. Y. Cai, Q. Ai, H. T. Quan and C. P. Sun: Sensitive chemical compass assisted 
by quantum criticality. Phys. Rev. A, 85 , 022315 (2012) 

[127] J. M. Cai and M. B. Plenio: Chemical compass model for avian magnetoreception 
as a quantum coherent device. Phys. Rev. Lett., Ill, 230503 (2013) 

[128] O. Efimova and P. J. Hore: Evaluation of nuclear quadrupole interactions as a source 
of magnetic anisotropy in the radical pair model of the avian magnetic compass. 
Mol. Phys., 107 , 665 (2009) 

[129] B. M. Xu, J. Zou, J. G. Li and B. Shao: Estimating the hyperfine coupling parame¬ 
ters of the avian compass by comprehensively considering the available experimental 
results. Phys. Rev. E, 88, 032703 (2013) 

[130] H. Lee, N. Yang and A. E. Cohen: Mapping nanomagnetic fields using a radical 
pair reaction. Nano Lett., 11 , 5367 (2011) 

[131] G. G. Guerreschi, M. Ticrsch, U. E. Steiner and H. J. Briegel: Optical switching 
of radical pair conformation enhance magnetic sensitivity. Chem. Phys. Lett., 572 , 
106 (2013) 

[132] G. Kodis, P. A. Moore, T. A. Moore and D. Gust: Synthesis and photochemistry of 
a carotene-porphyrin-fullerene model photosynthetic reaction center. J. Phys. Org. 
Chem, 17 , 724 (2004) 

[133] J. Cai: Quantum probe and design for a chemical compass with magnetic nanos¬ 
tructures. Phys. Rev. Lett., 106,100501 (2011) 

[134] A. E. Cohen: Nanomagnetic control of intersystem crossing. J. Phys. Chem. A, 113 , 
11084 (2009) 


56 


[135] A. W. Chin, A. Datta, F. Caruso, S. F. Huelga and M. B. Plenio: Noise-assisted 
energy transfer in quantum networks and light-harvesting complexes. New J. Phys., 
12 , 065002 (2010) 

[136] S. Shim, P. Rebentrost, S. Valleau and A. Aspuru-Guzik: Atomistic study of long- 
lived quantum coherences in the Fenna-Matthews-Olson complex. Biophys. J., 102 , 
649 (2012) 

[137] E. Rieper, Quantum Coherence in Biological Systems, Thesis, University of Braun¬ 
schweig, (2011) 

[138] J. C. Weaver and M. Bier, Bioengineering and Biophysical Aspects of Electromag¬ 
netic Fields-Chapter 7, Edited by F. S. Barnes and B. Greenebaum, CRC PRESS 
(2006) 

[139] J. C. Weaver, T. E. Vaughan and R. D. Astumian: Biological sensing of small field 
differences by magnetically sensitive chemical reactions. Nature, 405 , 707 (2000) 

[140] M. Tiersch and H. J. Briegel: Decoherence in the chemical compass: the role of 
decoherence for avian magnetoreception. Phil. Trans. R. Soc. A, 370 , 4517 (2012) 

[141] E. M. Gauger, E. Rieper, J. J. L. Morton, S. C. Benjamin and V. Vedral: Sustained 
quantum coherence and entanglement in the avian compass. Phys. Rev. Lett., 106 , 
040503 (2011) 

[142] P. Rebentrost, M. Mohseni and A. Aspuru-Guzik: Role of quantum coherence and 
environmental fluctuations in chromophoric energy transport. J. Phys. Chem. B, 
113 , 9942 (2009) 


57 


[143] J. Zhu, S. Kais, P. Rebentrost and A. Aspuru-Guzik: Modified scaled hierarchical 
equation of motion approach for the study of quantum coherence in photosynthetic 
complexes. J. Phys. Chem. B, 115 , 1531 (2011) 

[144] A. W. Chin, J. Prior, R. Rosenbach, F. Caycedo-Soler, S. F. Huelga and M. B. 
Plenio: Vibrational structures and long-lived electronic coherences. Nature Phys., 
9 113 (2013) 

[145] A. W. Chin, S. F. Huelga and M. B. Plenio: Coherence and decoherence in biological 
system: principles of noise-assited transport and the origin of long-lived coherences. 
Phil. Trans. Act. Roy. Soc., 370 , 3638 (2012) 

[146] F. Caruso, S. Montangero, T. Calarco, S. F. Huelga and M. B. Plenio: Coherent 
open-loop optiml control of light-harvesting dynamics. Phys. Rev. A, 85 , 042331 
( 2012 ) 

[147] V. Vedral, M. B. Plenio, M. A. Rippin and P. L. Knight: Quantifying entanglement. 
Phys. Rev. Lett., 78 , 2275 (1997) 

[148] S. Hoyer, A. Ishizaki and K. B. Whaley: Spatial propagation of excitonic coherence 
enables ratcheted energy transfer. Phys. Rev. E, 86, 041911 (2012) 

[149] R.Horodecki, P. Horodecki, M. Horodecki, and K. Horodecki: Quantum entangle¬ 
ment. Rev. Mod. Phys., 81 , 865 (2009) 

[150] H. J. Hogben, T. Biskup and P. J. Hore: Entanglement and Sources of magnetic 
anisotropy in radical pair-based avian magnetoreceptors. Phys. Rev. Lett., 109 , 
220501 (2012) 

[151] A. Peres: Separability criterion for density matrices. Phys. Rev. Lett., 77 , 1413 
(1996) 


58 


[152] M. Navascues, M. Owari and M. B. Plenio: Complete Criterion for Separability 
Detection. Phys. Rev. Lett., 103 , 160404 (2009) 

[153] G. Vidal and R. R. Werner: Computable measure of entanglement. Phys. Rev. A, 
65 , 032314 (2002) 

[154] S. Kais: Entanglement, electron correlation, and density matrices. Adv. Chem. 
Phys., 134 , 493 (2007) 

[155] K. Berrada, M. El Baz, H. Eleuch and Y. Hassouni: A comparative study of nega¬ 
tivity and concurrence based on spin coherent states. Int. J. Mod. Phys. C, 21 , 291 
( 2010 ) 

[156] Y. A. Lee and G. Vidal: Entanglement negativity and topological order. Phys. Rev. 
A, 88, 042318 (2013) 

[157] F. Verstraete, K. Audenaert, J. Dehaene and B. De Moor: A comparison of the 
entanglement measures negativity and concurrence. J. Phys. A: Math. Gen., 34 , 
10327 (2001) 

[158] X. Wang and S.-J. Gu: Negativity, entanglement witness and quantum phase tran¬ 
sition in spin-1 Heisenberg chains. J. Phys. A: Math. Theor., 40 , 10759 (2007) 

[159] M. Tiersch, G. G. Guerreschi, J. Clausen and H. J. Briegel: Approaches to measur¬ 
ing entanglement in chemical magnetometers. J. Phys. Chem. A, 118 , 13 (2014) 

[160] J. A. Pauls, Y. Zhang, G. P. Berman and S. Kais: Quantum coherence and entan¬ 
glement in the avian compass. Phys. Rev. E, 87 , 062704 (2013) 

[161] Y. Zhang, G. P. Berman and S. Kais: Sensitivity and entanglement in the avian 
chemical compass. Phys. Rev. E, 90 , 042707 (2014) 


59 


[162] R. Haberkorn: Density matrix description of spin-selective radical pair reactions. 
Mol. Phys., 32 , 1491 (1976) 

[163] N. Lambert, Y.N. Chen, Y.C. Cheng, C.M. Li, G.Y. Chen and F. Nori: Functional 
quantum biology in photosynthesis and magnetoreception. arXiv: 1205.0883 (2012) 

[164] R. Kubo: Stochastic Liouville Equations. J. Math. Phys., 4, 174 (1963) 

[165] A. T. Dcllis and I. K. Kominis: The quantum Zeno effect immunizes the avian com¬ 
pass against the deleterious effects of exchange and dipolar interactions. Biosystems, 
107 , 153 (2012) 

[166] 1. K. Kominis: Quantum Zeno effect explains magnetic-sensitive radical-ion-pair 
reactions. Phys. Rev. E, 80 , 056115 (2009) 

[167] G. E. Katsoprinakis, A. T. Dcllis and 1. K. Kominis: Coherent triplet excitation 
suppresses the heading error of the avian compass. New J. Phys., 12 , 085016 (2010) 

[168] 1. K. Kominis: Reactant-product quantum coherence in electron transfer reactions. 
Phys. Rev. E, 86, 026111 (2012) 

[169] J.N. Bandyopadhyay, T. Paterek, and D. Kaszlikowski: Quantum Coherence and 
Sensitivity of Avian Magnetoreception. Phys. Rev. Lett. 109 , 110502 (2012) 

[170] C. T. Rodgers, S. A. Norman, K. B. Henbest, C. R. Tinnnel and P. J. Hore: Deter¬ 
mination of radical re-encounter probability distributions from magnetic field effects 
on reaction yields. J. Am. Chem. Soc., 129 , 6746 (2007) 

[171] J. Cai, G. G. Guerreschi and H. J. Briegel: Quantum Control and Entanglement in 
a Chemical Compass. Phys. Rev. Lett., 104 , 220502 (2010) 


60 


[172] S. Maus, S. Macmillan, S. McLean, B. Hamilton, A. Thomson, M. Nair and C. 
Rollins, The US/UK World Magnetic Model for 2010-2015, NOAA Technical Report 
NESDIS/NGDC, (2010) 

[173] R. Mnlieim, Photobiology: the science of life and light. Lars Olof Bjorn (ed) 
Springer, New York, pp465-478 (2008) 

[174] Z. Huang and S. Kais: Entanglement evolution of one-dimensional spin systems in 
external magnetic fields. Phys. Rev. A, 73, 022339 (2006) 

[175] QC. Shi and S. Kais: Discontinuity of Shannon information entropy for two-electron 
atoms. Chem. Phys., 309, 127 (2005) 

[176] Z. Huang and S. Kais: Entanglement as measure of electronelectron correlation in 
quantum chemistry calculations. Chem. Phys. Lett. 413, 1 (2005) 

[177] M. Reed and B. Simon, Methods in Modern Mathematical Physics, Vol. 1, Sec. VI. 
6, Academic Press (1972) 

[178] G. Sadiek, Z. Huang, O. Aldossary and S. Kais: Nuclear-induced time evolution of 
entanglement of two-electron spins in anisotropically coupled quantum dot. Molec¬ 
ular Physics, 106, 1777 (2008) 


61 


